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Introduction 



■ Helicopter noise prediction is increasingly important 

>► certification 
>► detection 

■ A great deal of progress has been made since the mid 1980’s 

■ Purpose of this talk 

>- Put into perspective the recent progress 
>► Outline current prediction capabilities 
>► Forecast direction of future prediction research 
>► Identify rotorcraft noise prediction needs 
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Outline of Talk 



■ Introduction and Historical perspective 

■ Description of governing equations 

■ Current status of source noise prediction 

■ Future directions 

■ Summary 
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Rotor Source Noise 
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History of Helicopter Noise Prediction 


Propeller noise theory developed 
(steady loading, thickness) 


Importance of unsteady loading 
recognized 


Rotor noise theory development 
Helicopter rotor noise code development 

- excellent validation data available 

- large increase in computation power 


- 1940 - 
- 1950 - 
- 1960 - 
- 1970 - 
- 1980 - 
- 1990 - 
- 2000 - 


Helicopter rotor noise mechanisms 
proposed 


Ffowcs Williams-Hawkings equation 

- computer power limited 

- inadequate blade loading available 

(NR ) 2 program 

Kirchhoff formulation / 
quadrupole noise prediction / 
new application of FW-H equation 
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Available Methods for Rotor Noise Prediction 



■ Acoustic Analogy 

>► treats real flow effects by fictitious sources; exact in principle 
>► for rotor blades: Ffowcs Williams-Hawkings equation (1969) 

>- most developed, widely used in the helicopter industry 

■ Kirchhoff Formula 

>- originally suggested by Hawkings (1979); (Farassat and Myers 1988) 

>► method currently under development (development has been very rapid) 
>► depends upon high resolution aerodynamics input data from CFD. 

■ CFD based Computational Aeroacoustics (CAA) 

>- least mature 

> most computationally demanding 
>- advances in CAA will help other methods 
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■ Treats real flow effects by fictitious sources 

■ A mathematical device which is exact in principle 

■ Capable of supplying good qualitative and quantitative results 

■ For rotating blades 

>- Aerodynamic and acoustic problems separated 
>- Powerful methods of linear analysis can be used 
>► Inclusion of nonlinear effects feasible now 

■ Acoustic analogy is and will remain a very useful tool in 
aeroacoustics 
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Idea: rearrange governing equation into a wave equation 

continuity 


d_ 

dt 


dp + dpt. i ± = Q 


dt dx. 


d \ d pu i d _ 

d^\^r + d7, (pu ‘ Uj+p ^ =0 


d 2 p d 7 


dt 2 dxdx 

1 J 


(pu iUj +P iJ ) 


momentum (N-S) 


form wave equation 



where 


T u =pu i u j +P ij - c o p8 
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Ffowcs Williams-Hawkings Equation 
Derivation Procedure 



■ Embed exterior flow problem in 
unbounded space 

> define generalized functions 
valid throughout entire space 

> interpret derivatives as 
generalized differentiation 



r p 

/> 0 

p = < 

J Po 

/< 0 


pu { 

/> 0 

p«, = ■ 

0 

/<0 

p. = 

r p 

/>0 

V 

0 

V. 

/< 0 


■ Generalized conservation 
equations: 


dp 

dt 


+ 




+ P U i 



continuity 


dpu, t dpufij t &>> 

dt dxj dxj 


= ( P«i + (PUiUj + 


momentum 
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■ Numerical solution of the FW-H equation 

□ V(*. t) = |[P„ v„ 8(f)] - ^-[/,5(/)] + 

at ox. ox -ox. 

i i j 

■ Three source terms 

>- thickness source (monopole) 

- requires blade geometry and kinematics 

> loading source (dipole) 

- requires blade geometry, kinematics, and surface loading 

> quadrupole source 

- requires flow field (i.e., T { . ) around the blade (volume integration) 

■ WOPWOP+ implements all three of these source terms 
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Kirchhoff Derivation Procedure 


■ Use embedding procedure on wave equation 

>► define generalized pressure perturbation: 


\p' f> o 

|o /< o 


> use generalized derivatives 

>- generalized wave equation is Kirchhoff governing equation: 


I-.2 //■ — ,\ dp M dp d( ,M n \ d , 

□ p (x,t) = - — - + — p —0(f) - — (p n i o(f)) 

( at c on dt\ c ) dx ; 


V 

Qkir 
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Fgmiylgligf^gvglgBmgnl 


Model wave equation to solve (valid in entire unbounded space) 
U^(x,t) = Q(x,t)d(f) 


Integral representation of solution (Green’ s function 


f oo 


4 7I(l>(x, t) = | | 


Q{y,T)8(f)8(g) 


<5(g) 

4 nr 


) 


dydz 


Three potential formulations: 


4 = o 
8 = 0 


4 7t<j)(x,t)= f f cdTdt = f 1 

J J rsin0 J 


reGU) i 

dL = 

ret 

f 


A 

J 

/=o 

rl — M r \ 




-I ref 


collapsing sphere 
formulation 


emission surface 
formulation 


retarded time 
formulation 
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■ Retarded-time solution to FW-H equation (neglecting quadrupole) 



where Q = pv n and L = P ij h j 

■ Take derivatives inside integrals analytically (formulation 1A) 


4 7lp\x,t)- J 


/= 0 


Q + LJc 
r(l-M r ) 2 


+ 


I 

f = o 


dS+ J 

/= o 
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dS 
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NASA Rotor Noise Prediction Codes 



■ WOPWOP 

>► Uses FW-H equation, Farassat’s formulation 1A 
>► Used for discrete-frequency noise prediction 

> Representative of time-domain prediction codes (Primary U. S. code) 
>- Code features 

- Near and far-field acoustics 

- Forward flight and hover 

- Stationary and moving observers 

- Unsteady and impulsive loading allowed as input 

- Loading input may be analytical, computational, or experimental 

- Transportable, efficient, and robust 

■ WOPWOP+ 

>► includes a far-field quadrupole computation 
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NASA Rotor Noise Prediction Codes 



■ RKIR 

>► original code from Purdue University; modified by Sikorsky and 
NASA Langley to include all WOP WOP blade motions 

> utilizes Farassat and Myers’ Kirchhoff formulation for moving 
surfaces 

d p 

> require p, — ^ and Vp on the Kirchhoff surface 

ot 

■ FW-H/RKIR (prototype code) 

>► based on RKIR (Rotating Kirchhoff code) 

>► utilizes Farassat ’s formulation 1A (FW 
>- quadrupole source neglected; could be included 

■ Tiltrotor Aeroacoustic Codes (TRAC) 

>- collection of codes to predict the airloads, flow-field, and noise 
>► utilizes any of these codes to predict rotor noise 
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■ Thickness and Loading Noise 

■ Blade Vortex Interaction Noise 

■ High-Speed Impulsive Noise 

■ New Prediction Tools 

>► Kirchhoff Predictions 
>- FW-H Equation applied off the body (i.e. like a Kirchhoff formula) 

■ Broadband Noise 
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Predictions accurately reflect design changes 


llOkts 

upstream mic in TPP 
on advancing side 


ref: Brentner 1987 


Acoustic 




a) Rectangular planform 
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b) Tapered planform 
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Thigkngg^n^Lgg^in^Ngigg 



■ Predictions distinguish between small differences in input parameters 

■ Computations are efficient (29 CPU sec/observer on 22 M FLOPS workstation) 



ref: Brentner et al. 1994 



effect of pitching 
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effect of flapping 
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BVI Noise Prediction: with measuredairloads 


■ Amplitude, waveform, and spectra predicted well 

■ High temporal and spatial resolution of blade loads essential 




Harmonic number 


* microphone located upstream of rotor on advancing side, 25 deg. below TPP 

H = 0.152, C T /c 7 = 0.07, decent condition Ref: Brentner et al. 1 994, Visintainer et al. 1 993 
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BVI Noise Prediction: with calmtfatedairloads 


■ Near first principles prediction r e f: Tadghighi et ai. 1990 

■ Representative of state-of-the-art 


Comprehensive Code 


1 


CFD Flow Solver 


l 


Acoustic Prediction 
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■ High-speed impulsive (HSI) noise 

> particularly intense and annoying 

> occurs in high-speed forward flight 

> onset usually very rapid 

> primarily in-plane directivity 

■ HSI noise prediction 

> requires knowledge of 3D, nonlinear flow field 

> computationally intensive 

> modeled by FW-H quadrupole source 
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■ Importance of quadrupole source recognized 
Yu, Caradonna, and Schmitz (1978) 

>► simplified source strength 

>- far-field assumption / preintegration in z direction 
>- relatively immature flow field calculation 

■ Recent efforts 

>► Prieur (1986) - frequency domain, hover only 
>► Schultz and Splettstoesser (1987) - followed Yu et al. 

>► Farassat (1987-1991) - shock noise theory 

> Schultz et al (1994) - approx, source strength, both volume integration and 
preintegration 

>► Ianniello and De Bernardis (1994) - full volume integration 

■ NO readily available quadrupole prediction code in U.S. 
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■ Prediction by approximate quadrupole calculation 

>► Measured blade pressures and computed flow field used in prediction 


M h = 0.9 

hovering rotor 200 

mic in TPP 

Acoustic q 
pressure 

-200 

-400 

-600 


ref: Schultz and Splettstoesser 1987 
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Ref: Baeder 1991 


■ Prediction by direct CFD computation 

> Nonlifting, symmetric rotor in hover 


M h = 0.92 
hover 
mic in TPP 



Time, msec 
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■ FW-H source contributions linearly superimpose 

p\x,t ) = p'(x,t) + p' £ (x,t) + Pq (x,t) 

> develop quadrupole source prediction independently 

> can identify contributions from each source 

■ Current prediction codes based on FW-H equation 

> significant knowledge base 

> thickness & loading noise predictions very efficient 

■ Less demanding CFD computation 

> only compute the source region 

> don’t need to capture long-distance wave propagation 

■ Easy to study role of complicated rotor kinematics 
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l l J 

■ Source terms linearly superimpose 

p'(X,t) = p' t (X,t) + Pg(X,t) + Pg(X,t) 

■ Quadrupole source region is a volume 

>- needs large amount of data - 3D time dependent 
>- naturally separate 

■ Current WOP WOP very efficient 

>► desirable to not change thickness and loading now 

>► want to benefit from knowledge gained in thickness and loading noise 
development 


FW-H: 


□ 


p'(x,o = l-[p o v n 5(f)]- 


dx . 


+ 


d 


dx . dx . L l J 


T H(f) 
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Equation 

4 ltp' Q (x,t) = 


1 d 


-)2 ' 

bil 


T 


=/> o 


^dQdz 

r 


+ 


c dt 

d r r 37 - 7 


dt 


s 


rr u 


dCldz 


+c 


— °o />0 

371 - 7, 


n 

— oo f> 0 


rr u 


dQdz 



■ Interpretation 

>► f>0 - everywhere outside of blade surface 
>► dL 1 - element of collapsing sphere surface 

T,j = p Ui Uj +(p' - p'c 1 )5 ij 
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■ Assumptions 

>- Far-field observer 
>- In-plane observer 

■ Define new tensor 

oo 

Qu = J T iA 



■ Collapsing sphere approximated as a cylinder 

■ Integration in z is independent of source time 
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Fgi^Figl^AggigximgJjgn 



Contours of 
quadrupole 
source strength 
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WOPWQP+ Validation 



■ Validation case 

>► UH-1H, 1/7 th scale model rotor (untwisted) 

>► Experimental data available - Boxwell et al., Purcell 
>- Unique Euler calculation available (Baeder) 

- good resolution of flow field around blade 

- solution extends to microphone position at 3.09 R 

- symmetric solution 

■ Operating conditions for comparison 

> hover 

> M h ={0.88, 0.925} 

> inplane microphone at 3.09 R 
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Acoustic pressure, Pa 


■ Observer inplane at 3.09 R from rotor hub, M H = .85 




■ Quadrupole contribution roughly one-third that of thickness and loading 

■ Good agreement with Euler calculation and experiment 
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r Acoustic pressure, Pa 





■ Observer inplane at 3.09 R from rotor hub, M H = .88 




■ Good agreement with Euler calculation and experiment 


NASA Langley Research Center, Hampton, VA 


K. S. Brentner 


33 


■ Observer inplane at 3.09 R from rotor hub, M H = .90 




■ Quadrupole contribution is larger than thickness and loading and has 
steepened 

■ Retarded-time formulation does not allow all contributing panels to be 
included 
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■ Observer inplane at 3.09 R from rotor hub, M H = .925 




■ Quadrupole negative peak pressure shifts at higher speed 

■ Quadrupole contribution nearly twice that of thickness and loading 
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■ Contours of Q u 


M h =.88 



M h =.925 


UH-1 H model rotor 
untwisted; no lift 
hover test 




sonic circle 

extent of 

guadrupole 

integration 


WOPWOP+ preprocessor 
output 
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Acoustic pressure, Pa 


UH-1H Model Rotor Noise 

■ Components of acoustic pressure 




Time, msec 




Time, msec Time, msec 


UH-1 H model rotor 
untwisted; no lift 
hover test 




Thick + Load 

Quadrupole 

Overall 


Time, msec Time, msec 
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■ Observer inplane at 3.09 R from rotor hub, M H = .95 




■ Quadrupole term dominates pressure time history 

■ Predicted signal amplitude overpredicted 

■ Complete signal widening not predicted, but shock-like feature captured 
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Effigigng^ 



■ Preprocessor 

>► nominal run time: 3-5 CPU seconds 

■ Acoustic calculation 

> thickness and loading noise: - 5 CPU seconds 
>- quadrupole noise: - 11-17 seconds* 

> total: - 16-22 CPU seconds 

~ 45 CPU seconds when code forced to use 20pts/panel on last two rows 
CPU times for HP 735-99 scientific workstation 

■ Efficiency considerations 

>- quadrupole noise computation comparable to thickness and loading 
on a per panel basis 

>► adaptive quadrature enables use of a large number of quadrature 
points when needed 

>► reductions in CPU time possible 
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■ FW-H applied off the blade surface (like a Kirchhoff method) 

■ Kirchhoff method for moving surfaces 
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Not necessary to assume integration surface f=0 is coincident with 
body 


n 2 p\x,t) = 


dxdx 

1 J 

d 




dX; 


[{P lj n j +pii l {(ii ll - v„))3(/)] 


dt 
dx ; 


= — v. 


n i 


+ J t [( P° V n + P((«» - V „ )) 5 (/)] 


FW-H can be used as a Kirchhoff formula 


NASA Langley Research Center, Hampton, VA 


K. S. Brentner 


41 





■ Compare components from 
FW-H/RKIR with WOPWOP+ 

>► UH-1H rotor in hover 

> Hover solution from TURNS 
(Baeder) 

■ Two predictions necessary with 
FW-H/RKIR 

> thickness and loading from 
surface coincident with rotor 
blade 

> total signal (including 
quadrupole) from a surface 
approximately 1.5 chords away 
from blade. 

■ New application of FW-H 
equation retains advantage of 
predicting noise components 
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Manipulate FW-H source terms into form of Kirchhoff source 
terms (inviscid fluid) 


n 2 p\x,t) = Q kir + 


dxdx 

1 J 


[T,>W)\ 


d [pu i u j ]n,S(f) - f [pw,M„<5(/)] 


dx 


dx , 




■ Extra source terms are 2nd order in perturbations quantities 

■ FW-H and Kirchhoff source terms 

> equivalent in linear region ( p ' ~ c 2 p' u t « 1) 

> NOT equivalent in nonlinear flow region 
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■ UH-1H rotor 

> l/7th scale model 
>- untwisted blade 

■ Test setup (Purcell) 

> Hover, M H = 0.88 

> inplane microphone, 3.09 R 
from hub 

> minimal rotor lift 

■ Flow-field computation 

> full potential flow solver used 
(FPRBVI) 

> 80 x 36 x 24 grid (somewhat 
coarse) 

>- no rotor lift 
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Numerical Comparison: 


■ Principal advantage of the FW-H approach is insensitivity to 
surface placement 



time, msec 




data 


Kirchhoff 


(Note difference in pressure scales) 


NASA Langley Research Center, Hampton, VA 


FW-H 


K. S. Brentner 


45 





■ Advancing-side acoustic pressure ■ All three codes agree with each 

underpredicted other — non-lifting rotor 

■ Agreement with data is good 
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FW-H vs. Kirchhoff 



■ FW-H method of choice for aeroacoustic problems 

>► conservation of mass and momentum built in 

>► unified theory with thickness, loading, and quadrupole source terms 

> insensitive to integration surface placement 

■ FW-H approach the “better” than linear Kirchhoff because: 

> valid in linear and nonlinear flow regions 

>► surface terms include quadrupole contribution enclosed 
>► physical noise components can be identified with two surfaces 

■ The Kirchhoff approach 

> valid only in the linear flow region (not known a priori) 

- input data must satisfy the wave equation 

- wakes and potential flow field can cause major problems 
>► solution can be sensitive to placement of Kirchhoff surface 
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Broadband Noise 



■ Understanding 

>► Subjectively very important 

>► Many different mechanisms responsible - separate treatment for each 
>- Physical generation mechanisms well understood 

■ Prediction status 

>- Unsteady blade loads calculation difficult - classical methods used 
>► Frequency domain methods only - turbulence data in frequency domain 
>► Good prediction where turbulence statistics are known 
>► Good prediction of self-noise with semi-empirical methods 

■ Little explored approaches 

> Application of FW-H equation 

> Direct simulation of blade turbulence 
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Future Directions 



■ Ffowcs Williams - Hawkings equation 

>► Maturity level high — first choice for discrete frequency noise 
>► Efficient and robust codes currently available 
>- Solutions to current challenges in hand(BVI and HSI noise) 

■ Alternate approaches — feasible due to advances in CFD and 
computer technology 

> FW-H equation used as Kirchhoff method 
>► Direct computation of acoustics 

■ Relative importance of broadband noise increasing 

■ Continued work needed 

>- wake prediction 
>- aeroelastic coupling 

> full configuration aerodynamics/aeroacoustics 
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■ Rotor noise prediction capability is advanced 

>► Discrete frequency noise 

- Thickness and loading noise - prediction now routine 

- Blade- vortex interaction noise - good agreement demonstrated 

- High-speed impulsive noise - robust solutions available; 
depends upon CFD 

> Broadband noise 

- Semi-empirical predictions give good results for standard 
helicopter rotors 

■ Challenges for the future remain 

> Accurate prediction of high resolution airloads 

> Increased importance of broadband-noise prediction 

>► Systems noise prediction - component interaction; scattering; 
reflection 
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